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Data from ultrarelativistic heavy-ions collisions show evidence for temperature-fluctuations on 
the freeze-out surface of the expanding fireball. These may be remnants of density inhomogeneities 
in the initial collision overlap region. We present a power-spectrum analysis for heavy-ion collisions 
analogous to the analysis of the cosmic microwave background radiation. We use a Glauber model 
for eccentricity to extract the transfer-function needed to produce the observed spectrum and discuss 
its relation to the mean-free-path of the matter created in the collisions. 

PACS numbers: 



A commonly quoted goal of the heavy-ion programs at 
Brookhaven Lab (BNL) and CERN, is to recreate con- 
ditions similar to those shortly after the Big Bang when 
the universe was filled with a quark-gluon plasma (QGP). 
QGP can be created by smashing heavy nuclei together 
at relativistic speeds in collisions called " little bangs" . In 
the past few decades many theoretical and experimental 
advances have been made in the study of heavy-ion colli- 
sions. Comparisons to the early universe, however, have 
been scarce [2 [5] . ^ n this paper we explore an analogy 
between heavy-ion collisions and Big Bang cosmology. 
Using the heavy-ion equivalent of the map of the cosmic 
microwave background radiation (CMB), we determine 
the power-spectrum for the little bangs and we estimate 
the transfer-function necessary to produce the spectrum 
from the initial conditions. 

The hot fireball created in little bangs rapidly expands 
and cools and when cold enough forms hadrons. Even- 
tually, the system spreads out enough that hadrons stop 
interacting. This is called the surface of last scattering 
or freeze-out. The particles then free-stream to the de- 
tectors. From measurements of the number, mass, and 
momentum of these particles, we must infer the prop- 
erties of the matter that emitted them and extract the 
essential physics of the QGP. One of the most important 
discoveries at the Relativistic Heavy Ion Collider (RHIC) 
at BNL is that the tiny speck of QGP matter produced 
in the little bangs behaves much like a liquid [3 . This 
finding is based on the observation that the spatial asym- 
metries in the initial overlap zone show up as asymme- 
tries in the momentum distributions of final state par- 
ticles. The observed anisotropy is typically represented 
by the second Fourier component (z^) of the azimuthal 
distribution of observed particles relative to the reaction- 
plane @]. i>2 most strongly reflects the almond shape of 
the initial nuclear overlap region for non central colli- 
sions. The magnitude of «2 can be described surprisingly 
well with ideal relativistic hydrodynamic models suggest- 
ing a liquid-like behavior with minimal viscosity. The 



success of these models seems to indicate that the mean- 
free-path of interactions for the systems constituents is 
significantly smaller than the size of the system. 

Experiments at RHIC have also discovered correlations 
between particles that extend over a broad range in the 
longitudinal direction but are narrow in the azimuthal 
(transverse) direction forming a ridge [5]. A number 
of different scenarios have been proposed to explain the 
ridge |6t413j. These include minimum-bias (soft) jets in 
Au+Au collisions [BJ, soft gluons radiated by hard par- 
tons traversing the overlap region [7J, beam-jets boosted 
by the radial expansion [5], viscous broadening [9], and 
flux-tube like structures boosted by the radial expansion 
[HJ Q2] . The extent of the correlation in the longitudinal 
direction requires by causality that it must be established 
very early in the collisions [12] . One of us (PS) pro- 
posed that the correlation structures may be understood 
in terms of fluctuations of higher Fourier components of 
v„, particularly yj (vf ), that arise from anisotropics in the 
initial energy density converted into momentum space 
during the expansion [T3]. It was subsequently shown 
with the NEXSPHERIO hydrodynamic model, that in- 
deed, lumpy initial conditions lead to structures similar 
to those observed in the two-particle correlation measure- 
ments [14] . Alver and Roland [15] used the RQMD model 
to show that lumpiness in the initial collision geometry 
can lead to a yj in the azimuthal particle production. 
Petersen et al. ;16] carried out a similar analysis using 
an event-by-event hydrodynamic model. 

An analogy between the expansion of heavy-ion colli- 
sions starting from a lumpy initial energy density and the 
expansion of the universe starting with quantum fluctu- 
ations stretched to cosmological sizes was first pointed 
out by Mishra et al. [3J- They also proposed that the 
RMS values of v n (y/ (v%)) could be measured in heavy 
ion collisions analogous to the power-spectrum extracted 
from the CMB. They didn't however make a connection 
between yj (v%) and the already existing two-particle cor- 
relation measurements. In this work we use transverse 
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momentum {pr) correlations published by the STAR col- 
laboration [17] to extract the power-spectrum for heavy- 
ion collisions. Since the pr spectra reflects the temper- 
ature, pt correlations are sensitive to local temperature 
fluctuations. These measurements are directly analogous 
therefore, to the maps of the CMB. We use a Monte 
Carlo Glauber model [TH] for initial eccentricities to ex- 
tract the transfer-function required to convert the initial 
coordinate-space anisotropy into the anisotropy seen in 
momentum-space. This analysis facilitates a more direct 
comparison between relativistic heavy-ion collisions and 
the early universe. 

Analogy with Big Bang Cosmology: Measurements of 
the CMB reveal temperature fluctuations corresponding 
to over- or under-densities present at the surface of last 
scattering at about 300,000 years after the Big Bang [T5] , 
These density fluctuations ultimately explain the struc- 
ture in our universe (Fig. [Tjleft). Just as quantum fluc- 
tuations stretched to cosmic sizes by inflation show up 
in the CMB, we expect fluctuations from the beginning 
of the little bangs to show up in heavy-ion data (Fig. [I] 
right). Measuring temperature-fluctuations in the CMB 
required precise measurements at more than two million 
points in the sky. Enough photons are detected at each 
point to reconstruct the black-body spectrum from which 
the temperature is determined. In a heavy-ion collision, 
a few thousand particles are created at most, so a similar 
map cannot be made for each collision. But whereas we 
only observe one universe, billions of collisions are created 
in the lab. By studying px correlation data (sensitive to 
local changes in the p^-spectra and thus the tempera- 
ture) accumulated from millions of these collisions, we 
can search for evidence of hotspots on the surface of last 
scattering. 

THE UNIVERSE RHIC 




FIG. 1: Schematic of the expansion of the universe after 
the Big Bang (left) and the expansion of a fireball after little 
bangs (right). The illustration is by Alexander Doig. 

Survival of Density Fluctuations and Various Scales: 
In Fig. [2] we show the temperature fluctuations calcu- 
lated from a heavy-ion event generator near the begin- 
ning of the expansion and 4 fm/c later. We deter- 
mined these temperature profiles in the transverse plane 
at mid-rapidity by translating the energy-density pro- 
files of Werner et al. |20j into temperature, using the 
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FIG. 2: Temperature profile in the transverse plane for mid- 
rapidity at proper time r = 0.6 (left) and 4.6 fm/c (right). 

parametrized lattice QCD results for energy-density vs 
temperature from |21| . The simulations indicate that col- 
lisions of Au nuclei (12 fm across), may contain hotspots 
of size I spot ~ 1-5 fm and that remnants of those hotspots 
persists during the collisions evolution. We also consider 
the lengths of the acoustic horizon H and the mean-free- 
path l m f p of the systems constituents to be important. 

The acoustic horizon defines how far mechanical in- 
formation can have propagated through the medium at 
time r: H(t) = J^ fo c s (r)dT, where t/ is the freeze-out 
time. This relates to the growth rate of l spo t- We de- 
termined H{t) from lattice data on the speed of sound 
(c s ) vs energy density [22 and a hydrodynamic model 
to specify the energy density vs r [53] . Fig. [3] shows the 
acoustic horizon for QCD matter. The phase-transition 
from QGP to hadron-gas can be seen as a flattening in 
the slope of H at r « 10 fm/c when H is about 5 fm. 
The acoustic horizon also dominates the time dependence 
of the sound that an observer inside the medium would 
hear (see [53] and [55]). The sound at freeze-out is com- 
posed of a superposition of different waves with different 
frequencies that can be determined from the two-particle 
momentum correlation data. The horizon defines when 
frequencies can be heard: Only after half a wavelength 
fits inside the horizon would that wavelength become 
" audible" . This is the same effect that leads to the lack 
of large scale fluctuations in the CMB. 

The fact that hydrodynamic models do a reasonable 
job of predicting the value of vi suggests that l m fp can 
be considered small compared to the size of the sys- 
tem. By examining the power-spectrum of heavy-ion 
collisions which includes information for all values of n 
(beyond just n = 2 or n = 3), we hope to better con- 
strain l m f p . As we increase n, we reduce the length 
scale probed. We only expect an efficient conversion 
of coordinate-space anisotropies into momentum space 
when l m fp <C 2n(R)/n where (R) is the average radius 
of the systems constituents |26| . 

Power- Spectrum and Transfer-Function: We deter- 
mine the power-spectrum from two-particle momentum 
correlations (related to (p t i — (pt)){pt2 — (Pt))) vs the 
relative azimuthal angles between the particles Acj> [IT] . 
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FIG. 3: The acoustic horizon for heavy ion collisions at RHIC. 

A narrow peak positioned around small angle separation 
is observed in the data. This tells us that if a parti- 
cle comes out with above average pt, then the nearby 
particles also tend to have large pt- This is consistent 
with expectations from hotspots on the surface of last 
scattering. The correlation of these fast particles sug- 
gests that they are born out of the same high-density, 
high-temperature lump. We will not attempt to decom- 
pose the correlation into different components, i.e. jets 
and resonances and background. The power-spectrum we 
extract from data should and does contain all these con- 
tributions. We argue that jets do not dominate the ob- 
served correlation because the correlations are too large 
in magnitude, too narrow in </>, and too broad in rj: there 
will likely be some contribution but it is suppressed by 
1/multiplicity) . As for resonances, if a hotspot is there, 
it will emit massive and/or high momentum particles. 
The decay of a hotspot can proceed through decays into 
resonances. The power-spectrum reflects all these con- 
tributions. Our interpretation of the correlations data 
in terms of hotspots is supported by several pieces of 
ancillary evidence. 1) ii„-fluctuations are close to what 
we expect from density fluctuations from several models 
and the two-particle correlations data also match what 
we expect from these models [TH - fTSl |2"T] . This gives us 
confidence that the correlations are dominated by the 
over- and under-densities at the start of the expansion 
phase. 2) An improved description of particle pj- spectra 
is obtained when temperature fluctuations are considered 

ED- 

The pt correlations vs relative angles between the 
emitted particles (A(f> and Ar/) where parametrized in the 
STAR paper |17j . To extract the power-spectrum, we use 
that parametrization with Arj = and Fourier-transform 
the correlation function versus Atfi. The coefficients 
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a n = - /(A0)cos(nA0)d(A0). (1) 
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vs harmonic n make up the power-spectrum. If we had 
used number correlations instead of px correlations a n » 
Vn [25] ■ The power-spectrum for little bangs is shown in 
Figure [| (left). 

Having extracted the power-spectrum from the az- 
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FIG. 4: Left: The power-spectrum versus harmonic number 
at central rapidity. Right: Participant eccentricity for arbi- 
trary harmonic n for either perfectly central Au+Au collisions 
or the 5% most central. The intrinsic deformation of the Au 
nucleus is included. 



imuthal correlations of observed particles, we can com- 
pare that to the azimuthal distribution of matter in the 
initial overlap zone. This is calculated using the par- 
ticipant eccentricity {s ntPart ) for all harmonics n as in 
Ref. [15] . The fluctuations in the initial geometry cause 
the major axis of the eccentricity to fluctuate away from 
the reaction-plane direction. s njPar t is the eccentricity 
calculated along the major-axis. Fig. [1] (right) shows 
( £ n, P art) f rom a Monte Carlo Glauber model [T5] for per- 
fectly central collision (b=0 fm) and for the 5% most 
central collisions. The large n — 2 term persists even 
for central collisions because we include the intrinsic de- 
formation of the Au nucleus in our Monte Carlo. The 
n = 1 term is small because the participants are re- 
centered so that (x) — (y) — 0. For 6 = collisions 
we note that the eccentricity is nearly independent of n 
for n > 2. This is because for symmetric collisions and 
point-like participants, £ par t depends only on the number 
of participants, independent of n [30j . In this case, if all 
harmonics were converted into momentum space equally 
well, the final correlation function tend to a Dirac delta 
function at A0 = 0. We expect however, that the con- 
version of higher harmonic eccentricity will be damped 
due to the existence of the length scale l m fp- The con- 
version will be efficient only when l m f p < 2n(R)/n. We 
can investigate the damping of the higher modes by plot- 
ting the transfer-function which is the ratio of the power- 
spectrum in Fig. [Z] (left) to (e part n ) in Fig. [Z] (right) . 

Fig.pphows the transfer-function a n / (£ part , n )- We use 
( £ part,n) f° r a centrality range that matches the data. 
Understanding the shape of the transfer-function implies 
that one can understand the shape of the correlations 
data. The transfer-function shows that as n increases, 
the efficiency for converting coordinate-space anisotropy 
into momentum space quickly drops off as expected from 
the condition that the transfer-function should go to zero 
when l m f p ~ 2t:(R) jn. We can make a crude estimate for 
lmfp based on the transfer-function. If we take (R) = 3 
fm for the average radial position of a participant, and 
n = 5 as the harmonic beyond which the conversion is 
inefficient, then we get l m pf ~ 3.5 fm. This estimate 
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FIG. 5: The transfer-function for central Au+Au collisions 
extracted from p T correlations for n > 1 . 



corresponds to a viscosity five times larger than esti- 
mates based on the centrality dependence of v 2 [21] ■ Our 
crude estimate, however, is geometry based, not account- 
ing for the various phases of the expansion. The authors 
of Ref. [32] estimate the viscosity from the longitudinal 
width of the near-side peak in px correlations, while our 
estimate of l m f P is based on the transverse width. The 
transfer-function should be compared to a more complete 
model of heavy-ion collisions in order to understand the 
effects of the expansion velocity, viscosity in the QGP, 
and viscosity in the hadronic phase. 

In this letter we considered the acoustics of heavy ion 
collisions and discussed the analogy with the early uni- 
verse. We presented the power-spectrum from heavy 
ion collisions and derived the transfer-function needed 
to convert spatial correlations from the initial conditions 
into the px correlations measured at RHIC. We find the 
transfer-function required to describe the RHIC data can 
be easily understood in terms of an inefficiency in conver- 
sion of coordinate space anisotropy into momentum space 
when l m fp > 2n(R)/n. We used this transfer-function 
to make a rough estimate of the mean-free-path of the 
systems constituents. This approach represents a new 
method for determining the characteristics of heavy-ion 
collisions and the QGP and should be further investi- 
gated. 

Acknowledgments: While preparing this draft, a pa- 
per on a similar topic was posted was posted [33]. We 
thank Marcus Bleicher and Klaus Werner for providing 
the data for the energy density profile and Alex Doig for 
preparing several illustrations. AM thanks FIAS for the 
kind hospitality. 



[2] A. P. Mishra et al., Phys. Rev. C77, 064902 (2008). 
[3] M. Gyulassy, L. McLerran, Nucl. Phys. A750, 30-63 
(2005); 

[4] The reaction-plane is defined by the beam axis and the 
vector connecting the centers of the colliding nuclei. 

[5] M. Daugherity [STAR Collaboration], J. Phys. G 35, 
104090 (2008); J. Adams et al. [STAR Collaboration] 
Phys. Rev. Lett. 95:152301, (2005); J. Adams et al. 
[STAR Collaboration], Phys. Rev. C 73 (2006) 064907; 
B. I. Abelev et al. [STAR Collaboration], Phys. Rev. C 
80, 064912 (2009); J. Adams et al. [Star Collaboration], 
Phys. Rev. C 75 (2007) 034901; B. Alver et al. [PHOBOS 
Collaboration], J. Phys. G 35, 104080 (2008). 

[6] T. A. Trainor [STAR Collaboration], AIP Conf. Proc. 
828, 238 (2006). 

[7] A. Majumder, B. Muller, S. A. Bass, Phys. Rev. Lett. 
99, 042301 (2007). 

[8] E. V. Shuryak, Phys. Rev. C76, 047901 (2007). 

[9] S. Gavin, G. Moschelli, J. Phys. G G35, 104084 (2008). 
[10] S. A. Voloshin, Phys. Rev. B 632, 490 (2006). Nucl. Phys. 

A 802, 107 (2008). 
[11] S. A. Voloshin, Nucl. Phys. A749, 287-290 (2005). 
[12] A. Dumitru et a/.,Nucl. Phys. A810, 91 (2008); S. Gavin, 
L. McLerran, G. Moschelli, Phys. Rev. C79, 051902 

(2009) . 

[13] P. Sorensen, J. Phys. G: Nucl. Part. Phys. 37 094011 

(2010) . 

[14] J. Takahashi et al.,Phys. Rev. Lett. 103, 242301 (2009). 
[15] B. Alver, G. Roland, Phys. Rev. C81, 054905 (2010). 
[16] H. Petersen et al, arXi v:1008.0625| [nucl-th]] . 
[17] J. Adams et al. [STAR Collaboration], J. Phys. G 32, 
L37 (2006). 

[18] M. L. Miller et al, Ann. Rev. Nucl. Part. Sci. 57, 205 
(2007). 



[19 
[20 
[21 

[22 

[23 
[24 
[25 
[26 
[27 



[28 
[29 



[1] K. Yagi, T. Hatsuda and Y. Miake, Camb. Monogr. Part. 
Phys. Nucl. Phys. Cosmol. 23, 1 (2005); H. Schade and 
B. Kampfer, Phys. Rev. C 79, 044909 (2009); D. Boy- 
anovsky, H. J. de Vega, D. J. Schwarz, Ann. Rev. Nucl. 
Part. Sci. 56, 441-500 (2006). 



http://map.gsfc.nasa.gov/ 
K. Werner et al. |arXiv:1004.0805| [nucl-th]. 
A. Bazavov et al Phys. Rev. D 80, 014504 (2009); 
S. Borsanyi et al, |arXiv:1007.2580 ||j [hep-lat], 
P. Huovinen, P. Petreczky, Nucl. Phys. A837, 26-53 
(2010). 

P. F. Kolb, Heavy Ion Phys. 21, 2 43 (2004). 
http://soundofthelittlebang.com/ 
http://www.youtube.com/watch?v=jF8Q03Cou-Q 
n = 1 is suppressed by momentum conservation. 
P. Sorensen, Proceedings of 24th Winter Workshop on 
Nuclear Dynamics, South Padre Island, Texas, 5-12 Apr 
2008. 

Z. Tang et al, Phys. Rev. C 79, 051901 (R) (2009). 
T. A. Trainor and D. T. Kettler, In t. J. Mod. Phys. E 
17, 1219 (2008). T [arXiv:0704.1674| [hep-ph]]. 
[30] Giving the participants a size causes (£n, P art) to drop as 
a function of n. The drop occurs at n ~ tt{R) /r part where 
r P art is the radial size of the participant and (R) ~ 3 fm 
is the average radial position of a participant. This drop 
is similar to the drop we expect in the transfer-function 
due to lmfp and the effects could be confused in the data. 
[31] H. J. Drescher et al., Phys. Rev. C 76, 024905 (2007). 
[32] S. Gavin and M. Abdel-Aziz, Phys. Rev. Lett. 97, 162302 
(2006). 



[33] P. Staig and E. Shuryak, arXiv: 1008.3139 [nucl-th] 



